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By inelastic scattering of polarized neutrons near the (200)-Bragg reflection, the susceptibilities 
and linewidths of the spin waves and the longitudinal spin fluctuations, 5Ssm(q) and (5S3(q) |[ Ms, 
respectively, were determined separately. By aligning the momentum transfers q perpendicular to 
both SSsw and the spontaneous magnetization Ms, we explored the statics and dynamics of these 
modes with transverse polarizations with respect to q. In the dipolar critical regime, where the 
inverse correlation length Kz{T) and q are smaller than the dipolar wavenumber qd, we observe: (i) 
the static susceptibility of (5Sj„(q) displays the Goldstone divergence while for 5Sj(q) the Ornstein- 
Zernicke shape fits the data with a possible indication of a thermal (mass-) renormalization at 
the smallest q- values, i.e. we find indications for the predicted 1/q divergence of the longitudinal 
susceptibility; (ii) the spin wave dispersion as predicted by the Holstein-Primakoff theory revealing 
qd = 0.23(1) in good agreement with previous work in the paramagnetic and ferromagnetic 

regime of EuS; (iii) within experimental error, the (Lorentzian) linewidths of both modes turn out 
to be identical with respect to the g^-variation, the temperature independence and the absolute 
magnitude. Due to the linear dispersion of the spin waves they remain underdamped for q < qd- 
These central results differ significantly from the well known exchange dominated critical dynamics, 
but are quantitatively explained in terms of dynamical scaling and existing data for T > Tc- 
The available mode-mode coupling theory, which takes the dipolar interactions fully into account, 
describes the gross features of the linewidths but not all details of the T- and q-dependencies. 

PACS numbers: 68.35.Rh, 75.40.Gb 



I. INTRODUCTION 



Neutron scattering has been demonstrated as an ex- 
tremely useful probe of the spin fluctuations (SS((7l[aear 
the Curie temperature of Heisenberg ferromagneta^lia'&Q. 
Under the assumption that the isotropic exchange in- 
teraction dominates the ordering process, the early re- 
sults could be welUnterpreted in terms of the dynamical 
scaling hypothesiaa. In its simplest form this hypothe- 
sis states that as for the static spin correlations also the 
temperature and q-dependencies of their characteristic 
frequencies are described by homogeneous scaling func- 
tions that depend only on a single variable, q^{T), where 

denotes the correlation length of the order parameter 
fluctuations. 

Signatures of the inevitable, anisotropic dipole-dipole 
interaction on the fluctuations of Heisenberg ferromag- 
nets have been first detected by measurements of the 
relaxation rate r(0) of the homogeneaus (5S'(0) mode in 
the paramagnetic phase of CdCr2Se4 a and subsequently 
also, of Fe by neutron spin-echo measurements at small 
q Id. Based on these signatures and also on first theo- 
retical approaches, which treated the dipolar interaction 



as a perturbatioB-pf the isotropic fluctuationsErB, it has 
been conjecturedt2l that the dipolar forces should gain im- 
portance somewhere in the socalled dipolar critical (DC) 
regime, where = k and q arc small compared to the 
dipolar wavenumber (see Fig.^). For a given ferro- 
magnet, this quantity measures the strength of the dipo- 
lar interaction relative to the exchange interactions. It 
has been introduced in renormalization gmup (RG) cal- 
culations of the static correlation functionsU, recognizing 
that the dipolar anisotropy breaks the rotational invari- 
ance of the fluctuations (5S(q) with respect to q. The 
demagnetizing effect of the dipolar interaction on the lon- 
gitudinal modes 6S^ \\ q, prevents them from becoming 
critical, while the remaining two transverse modes (see 
Fig.|l|c) are driving the ferromagnetic transition. At first 
these dipolar anisotropic fluctuations have been realized 
by elastic scattering of polarized neutrons pabove Tc of 
the Heisenberg ferromagnets EuO and EuSIIj, where for 
the latter it was also possible to measure directly the 
characteristic wavenumber, qd — 0.22(5)A^^. 

On the theoretical side, the implications of the dipo- 
lar critical fluctuations on their dynamics have been 
fully taken into account only by the mode-mode coupling 
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FIG. 1: a) Dipolar regime for static critical behavior below 
Tc explored in the present work (shaded area), with qd — 
0.23 A"^ for EuS. The dotted lines define the exchange (EH) 
and dipolar (DH) hydrodynamic critical regimes above and 
below Tc, respectively, b) Spin fluctuation modes deflned 
by the geometry of our polarized neutron experiment below 
Tc of EuS, where the momentum transfer q = Q — T200 is 
kept perpendicular to the spontaneous magnetization. Spin 
wave (SSsw) and longitudinal {SSz) fluctuations are detected 
separately by the spin-flip and non spin-flip intensities for 
neutrons with incident polarization CTn and wave vector Q. c) 
Definition of the magnetic modes with respect to the reduced 
momentum transfer q 



(MMC), approachliJoliMX Above Tc, rather convincing 
agreement was obtainedllallil for the critical slowing down 
of F" (q) observed for g ^ and-^ Tc on the trans- 
verse fluctuations of FeQ and EjiSIIZI as well as for the lon- 
gitudinal ones {a = L) of EuSEZl. Additional strong sup- 
port for the MMC-results came from quantitative anal- 
yses of the relaxation rates r^(q, T) in the two limiting 
cases (J = 0, T > Tc and 31 = Tc, 9^0 for the archetype 
Heisenberg ferromagnetsH As one of the striking results 
we mention the crossover from F"^(g > qd,Tc) ~ q^^"^ in 
the exchange dominated-cefficne to r'^(g <C qd,Tc) ^ q^ , 
deep in the dipolar oncoEjO. The latter behavior cor- 
responds to the conventional (van Hove type) slowing 
down, characterized by a non-critical Onsager co efficien t 
of the spin dynamics, L-^(q, T) = F^(q) • X"^(q) li^'tllllS, 
where x^(q) is the static susceptibility of the transverse 
fluctuations. It was shown that this central quantity de- 
pends in a universal manner only on Tc and qd, to be 
discussed in Sect.^. 

Below Tc, the situation becomes more complicated 



because the symmetry of the fluctuations is further re- 
duced by the appearance of the order parameter, i.e. 
the spontaneous magnetization M^. As illustrated in 
Fig.^, it is common sense to distinguish there between 
the spin waves (5S<;tu(q) _L M^Gz and the longitudinal 
modes (5S2(q) || MsGz- A systematic classification of 
critical behaviors in the q — K(r)-plane, which considers 
the Heisenberg exchange and the dipolar interaction on 
an equal footLag, has recently been performed by Schinz 
and Schwabllla. According to their results, we have de- 
picted the dipolar critical (DC) regime in Figj^a, which 
apart from a dipolar hydodynamic (DH) regime at rather 
small q is bound from above by the dipolar wavenumber 
qd in q-direction as well as in the K2(T)-direction. 



Since among the ar,dj,el^pical Heisenberg ferromagnets 
experimental studies of dipolar effects^lill^l 



is largest for EuSll3E3, this material 



ered for 
3. Here 

we report results of a first systematic study of the mag- 
netization dynamics in the DC-region, which covers the 
shaded area in Fig.|l|a. Our principal goal is to determine 
the susceptibilities and linewidths of the spin waves and 
of the longitudinal fluctuations and to examine their q- 
and temperature variations. The experimental access by 
means of inelastic scattering of polarized neutrons around 
a finite Bragg peak is described in Section |l[ As illus- 
trated by Fig.^ we choose a configuration where the 
momentum transfer q occurs in directions perpendicu- 
lar to the order parameter Ms, which allows to define 
and determine the transverse polarizations of the spin 
wave and of the longitudinal modes with respect to q. 
This will turn out to be essential for the discussion by 
means of dynamical scaling, because the appearance of 
the order parameter does not lift the symmetry with re- 
spect to q so that both modes retain their transverse 
polarization from above Tc- In Section [11, we give some 
examples for constant Q-scans and the analysis of the 
inelastic cross section. We evaluate there the relevant 
observables, i.e. static susceptibilities, the spin wave fre- 
quencies and the linewidths. Their detailed temperature 
and q- variations are presented in Sect. [V. In particular. 



we determine there the dipolar wavenumber qd from spin 
wave energies and the correlation length £,{T) — k~^{T) 
of the longitudinal fluctuations from their static suscep- 
tibilities. Backed with these findings, we discuss the cen- 
tral results, i.e. the damping of the spin waves and of the 
longitudinal fluctuations, in Sect.^. Here the objectives 
are twofold. First, basing on the existing data for the 
paramagnetic side of EuS we want to examine whether 
and how the scaling hypothesis, which is rather special 
for the dipolar interactiond, works, and second we will 
compare our results to the rather detailed predictiaas 
of the recent numerical solutions of MMC equationdiS. 
To produce explicit values for the linewidth, this MMC 
approach had to introduce several assumptions which, of 
course, need to be checked by the experiment. The paper 
closes with a brief summary and outlook. 
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FIG. 2; Rocking curve of the T200 Bragg peak of the isotopi- 
caUy enriched sample ^^"^EuS that is composed of more than 
100 individual single crystals. The mosaic of the sample is 
77 = 0.75° ±0.02°. 



II. EXPERIMENT 

The neutron scattering experiments were performed 
using the triple axis spectrometer IN14 at the ILL in 
Grenoble, with polarization analysis. The isotopically 
enriched sample ^^'^EuS was composed of roughly 100 
single crystals, aligned such that the overall mosaic was 
7] ~ 0.75° (see Fig. ||). In EuS the magnetic i?-u^+ ions 
{S — 7/2,gL — 2) are arranged on a face centered cu- 
bic lattice, ao = 5.973A. The sample was mounted in- 
side a superconducting magnet providing vertical fields 
By up to 4.5 T. All measurements were conducted near 
the (200) Bragg peak with q along the [hOO] direction 
(see Fig.nb) using neutrons with fixed incident energies 
= 3.0 and 3.8 meV and coUimations 37'-37'-40' along 
the beam direction from the monochromator to the detec- 
tor. Due to the mosaic structure of the sample shown in 
Fig.^ the measurements in the other q-direction suffered 
from a high background from the elastic Bragg intensity. 

The neutron cross section for neutrons scattered from 
an isotropic magnetic material is given byE3 



do- kf 
dVLdE' ^ "fc" 
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iVF(Q)2exp(-2T4^(Q)) 
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)2^)5„Mq:^)- (1) 



Q = kj — k/ and Tilo — Ei — Ef are the total momen- 
tum and energy transfers from the neutron to the sam- 
ple, respectively, where Eij = {Tikf ^ / {2nnn) is the neu- 
tron energy, a and (3 designate the cartesian compo- 
nents x,y,z, (7ro/2)2 = 72.65 • 10"^ barn//i|, F is the 
magnetic form factor for the Eu^+ ion, the exponential 
term is the Debye- Waller factor, Q = Q/|Q|, and N 
is the number of the magnetic ions in the sample. The 
reduced momentum transfer is defined by q = Q — T200, 
where T200 is the position of the nearest Bragg peak. 



"''200 = (200)27r/ao in our case, (Fig. |^b). 

The a/3-component of the scattering function 
S'Q^(q, w) is related to the imaginary part of the 
susceptibility by the fluctuation-dissipation theorem 

'S'a/3(q,t^) = (n + l)-3Xa/3(q,'^), 
TT 

where the thermal population factor is given by (n) = 
[exp(?iw/fcBT) — 1]^^. Our experiments have been con- 
ducted close to the Curie temperature Tc = 16.25 K 
and at energy transfers Tilo <^ ksT , so that Sap becomes 
directly proportional to T'^Xap 



5*0/3 (q,w) 



ksT 1 

-T -^Xa/3(q,t^) 



(2) 



Thus, the neutron scattering cross section reflects di- 
rectly the q- and w-dependence of the susceptibility com- 
ponents of the sample. 

According to Eq.(|l]), the neutrons couple only to spin 
fluctuations 6S that are perpendicular to the scattering 
vector Q = T2oo+q- Therefore, the cross section contains 
contributions from the longitudinal fluctuations (5S^ and 
from the spin wave scattering that corresponds to ex- 
citations with (5S^„ perpendicular to M. Both modes 
have a transverse polarization with respect to q, and in 
the following we omit the superscript T. The fluctua- 
tions parallel to M can be separated from the spin wave 
modes by analyzing the polarization c/ = if .; of the 
scattered neutrons with respect to the incident polar- 
ization, (Ti II B„. As can be inferred from Figj^b, the 
parallel fluctuations give rise to non-spin flip scattering 
and the spin wave modes give rise to spin flip scattering. 
Therefore, the susceptibilities Xz(q, w) and Xsu)(q, of 
both magnetic modes, can be determined unambiguously 
in a vertical field Bt,. In order to minimize the influence 
of the magnetic field on the spin fluctuations as much as 
possible, we adjusted B^, so that the internal magnetic 
field B = B^ - N^M{T), with = 0.05, was smafl, 
however large enough to remove the domain walls. The 
flipping ratio was between 2.1 < R < 3.5. 

In a first step, all the measured data I^^^ was corrected 
for the finite fiipping ratio R according to: 
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where I^'^ designates the scattered intensity from the po- 
larization fj, to the polarization v of the incident and scat- 
tered neutrons, respectively. In a second step, a back- 
ground being determined in the paramagnetic phase at 
T = 80 K > Tc and in the ordered phase at T = 1.78 K, 
was subtracted from the data. The energy-independent 



4 




-0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 

Energy Transfer (meV) 

FIG. 3: Constant-Q scans for a) spin flip and b) non-spin 
flip neutrons scattered from spin wave and longitudinal fluc- 
tuations, respectively, at T = 0.92Tc = 15.0 K. c) Spectrum 
for unpolarized neutrons, calculated from a) and b), showing 
the depression of the longitudinal contribution. 

contribution was 1 count/8.3 min for.. all the measure- 
ments in the range 0.06 < ^ < 0.18 oEj, where the re- 
duced momentum transfer is measured in reciprocal lat- 
tice units C = g/(27r/ao). The peak intensity of the elas- 
tic background was 4 counts/8.3 min and 8 counts/min 
for the spin-flip and non-spin-flip scattering at the posi- 
tion T + q, respectively. 



III. CONSTANT-Q SPECTRA 

The inelastic magnetic scattering has been determined 
at several temperatures and momentum transfers by per- 
forming constant-Q scans. Fig. ^ shows the cross sections 
at C = 0.18 as measured along the [100] direction at the 
(2 0) Bragg reflection at T = 15 K in a field = 100 
mT. The spin-flip data clearly reveals spin waves, while 
the non-spin flip data is quasielastic and has a width 
(half-width at half maximum) that is roughly a factor 
of two smaller than the energy of the spin waves. Sim- 
ilar data was collected for many different temperatures 
1.78 < T < 80 K and in appropriate flelds 30 mT < By < 



500 mT, as described in Section ||. 

Scattering of polarized neutrons is a very efficient way 
to properly separate the longitudinal from the transverse 
excitations in an Heisenberg ferromagnet. This is demon- 
strated in Fig. 1^ where we show for comparison the 
spectrum calculated from the (flip efficiency) corrected 
intensities of Fig. ||a,b: 

2 1 

as it would be measured by means of unpolarized neu- 
trons. The comparison shows clearly that the longitu- 
dinal scattering cannot be reliably extracted by means 
of unpolarized neutrons because the spectral widths of 
Isw and Iz are similar. Moreover, reliable positions and 
widths of the spin waves can only be obtained if polarized 
beam data is used. 

In order to analyze the data we have employed the 
scattering function, Eq.(|l|), which via Eq.(||) is directly 
related to the dynamical susceptibilities Xz and Xsw- To 
allow for a cpmparison between our results and the ex- 
isting theoryll3, we assume Lorentzian spectral weight 
functions for the dynamic susceptibilities of the modes 
fi — ZjSW and obtain 

^^^^ \ 
^(c. + c.^(q))2-f (r,)^;' ^ > 

where hujsw = Eswitl) is the spin wave energy and r^(q) 
are the linewidths. As the scattering from the lon- 
gitudinal fluctuations proves to be quasielastic, we set 
ojzifl) = 0. The data was fitted by convoluting the scat- 
tering functions 5'^(q, w) with the four-dimensional reso- 
lution function of the spectrometer to give //^(q, uj). The 
three parameters iS^iii, T^, X;j(q) and a common normal- 
ization parameter were varied for each constant-Q scan 
such that x^ "^as minimized. The solid lines in Fig. ^,b 
are fits to the data using Eq. (^) . They describe the data 
well. Figure |^ shows that the width of the longitudinal 
fluctuations Tz increases with inajeasing q in qualitative 
agreement with existing theoriesElclEJ. A more detailed 
comparison will be performed in the discussion. 



IV. ANALYSIS OF THE SPECTRA 

From the fits of the spectra using the double Lorentzian 
scattering function, Eq. as described in the previ- 
ous paragraph, we extracted the static susceptibilities 
Xp(q,t^ = 0, T) = X/i(q,T'), the energies Esw{q) and 
linewidths HTswi^} of the spin waves, and the linewidth 
hTziq) of the quasielastic scattering. These quantities 
will be discussed in the following. 
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FIG. 4: Spectra of the longitudinal fluctuations recorded at 
T = 0.96Tc and different momentum transfers q. 
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A. Static Critical behavior 

In order to compare the experimentally determined 
static susceptibilities with the theory we refer to the 
expressipHS from the theoretical work by Schinz and 
Schwabij who have presented x^(q, T) for all q- values 
and temperatures below Tc- Especially, one finds for the 
susceptibilities Xsw and Xz, deep in the dipolar regime as 
covered by our experiment (Fig.|l|a): 



+ Ki(q,T) 



(4) 



(5) 



The q^^-divergence is characteristic of gapless spin 
waves with transverse polarization 5Ssw -L q, also re- 
ferred to as Goldstone modes, while the fluctuations par- 
allel to M acquire a (thermal) mass, k^_{T). At small 
momentum transfers g <C this term is renormalized 
by the spin wave fluctuations and becomes q-dependent 



FIG. 5: a) Inverse of the integrated intensities of the spin 
wave and longitudinal spin fluctuations versus , fltted to the 
inverse static susceptibilities, Xb-^ ^ il^ + '^9) Xz^^il ~^ 
0) ~ (g^+K?. (T, g)), Eq.(|^ b) Temperature dependence of the 
inverse correlation length of the longitudinal magnetization 
fluctuations K^(r) = C^{T), defined by Eq.(|a). 



' 18 1 



aKz{T)/q' 



(5a) 



where — denotea-l|ii« inverse correlation length 
below Tc and a ~ 2/9 ll§E3. One consequence of this 
effect has been emphasized for the homogeneous suscep- 
tibility alreadjtiii the original spin wave work by Holstein 
and PrimakoflcZl. They predicted the singularity for van- 
ishing magnetic field, Xz{<l = 0,T, B ^ 0) ^ B^^'^, which 
in fact has been confirmed in experiments on EuS and 
EuOEa. More recently, the crossover from the Ornstein- 
Zernicke type behavior of the zero-field susceptibility to 
^ Xzil Kz) q^^ singularity has been obtained by 
the RG theoryca, but a clear experimental evidence is 
yet lacking. 
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FIG. 6: Dispersion of the spin waves with transverse polar- 
ization (55^ (q) (see Fig.^). The energies have been nor- 
mahzed by the reduced spontaneopsi magnetization nis — 
Ms{T)/MsiO) = 1.18(1 - T/Tcf-^^ M a.wA fitted to the pre- 
dictions of the Holstein-Primakoff theoryE3,Eq.(^). 



A summary of our static results are depicted in Fig.|^a, 
where the inverse (w-integrated) intensities of the spin- 
flip and non-spin-flip channels measured at five differ- 
ent temperatures are plotted against q^. The former 
ones turn out to be independent of temperature and 
display a clearcut g^-dependence as predicted by Eq.(||) 
for the spin wave susceptibility Xsw- Note that by us- 
ing a representation of q^, which involves the spin wave 
stiffness Dsw{T)q^ = gLfJ-sMsiT), see Ref. |l^ and 
Eq.(7) below, one finds an equivalent form, Xsw {q ' — ^ 
0) = gLfiBMs{T)/D,UT)q^- In Fig.|a, at the lowest 
a slight offset is seen, Xsw{q = 0) — ^^/k^ with 
Kg — O.OAA^^. We ascribe it to the presence of a small 
gap in the spin wave spectrum, which may be associated 
with the small cubic anisotropy of EuS and the finite 
internal field required to remove the domains. 

In contrast to Xswi the susceptibility of the parallel 
fluctuations exhibits a strong temperature dependence. 
Figure ||a displays the fits of the inverse intensities to 
Eq.(^. At large q^, the Ornstein-Zernicke behavior is 
obeyed which by Eq.(5a) defines the inverse ferromag- 
netic correlation length Kz{T) = ■\/l8/29 k- (q » Kz,T), 
depicted in Fig.p3. Obviously, the temperature depen- 
dence can be well described by the critical law, 



K,{T)^K,{0){1-T/Tcr 



(6) 



with Tc = 16.25(5)K. The critical exponent i^' = 0.68(2) 
agrees wifh the value obtained above Tc of EuS, v — 
0.70(2) c3, as predicted by the static scaling hypothesis. 
For the critical amplitude we obtain Kz{0) = 0.91(5)A^^. 
Comparing this value to the amplitude of the param=. 
agnetic correlation length of EuS, Kp(0) — 0.53A~^ E£l 
we find Kz — l-7Kp, which is bracketed by the mean 
field value, \/2/pto) and 2.02k(0) obtained by consider- 
ing fluctuationsEj. 

At rather small q, we observe a downward bending 



of I~^{q), which we try to associate with the crossover 
X^^i<l ^ q following from Eq.(5a). Though 

the errors are fairly large, we fitted the inverse inten- 
sities to this prediction and found a=0.20(5), which 
agrees surnrisiugly well with a=2/9 predicted by various 
approacheslljcj. Regarding the fact that the tempera- 
ture variation of this (q)-hendmg is well reproduced, 
we believe that this constitutes the first signature of the 
^"^-singularity oixziq) induced by the Goldstone modes. 
Interestingly, thispwariation should also hold for the dipo- 
lar critical regimelij. 



B. Spin wave dispersion 

The spin wave energies, normalized to the reduced 
spontaneous magnetization nisiT) — AIs{T)/Ms{0) = 
1.18(1 - r/Tc)°-^^BK are shown in Fig.|. For the 
present orientation of q being perpendicular to M^., spin 
waveEH and linear responscEj theories predict for the q- 
dependence: 



E{q)^Doms{T)q^l + 



2\ 1/2 



(7) 



As for the static susceptibility we admitted the pres- 
ence of a small gap q"^ = q'^ + nj.^. Due to the well- 
known depolarization of one of the precessing compo- 
nents (5Ssu,, the dipolar wavenumber enters Eq.(^ via 
q^ = gi^fj,BMs{0)/ Dq to cause a crossover from the 
quadratic dispersion ai q ^ qd to the linear law for 
qd ^ q > The solid line in Fig.|| represents the 

fit of the data to Eq.(0) with ^ 3.02 meVA^.qd = 
0.23(1) and Ksw = 0.01 A"^. By comparing this 
fitted dipolar wavenumber to q^ = 0.22(5) A"jias deter- 
mined from paramagnetic neutron scattering)13 suggests 
that this quantity is not renormalized by critical fiuctu- 
ations when passing Tc- Sucii-|effect has been conjec- 
tured by Fisher and Aharonytil. However, some indi- 
cation for the absence of such renormalizations through 
critical fluctuations in the dipolar regime has alreajiy 
been realized during a previous determination of qj^. 
There, using the mean field expression for the critical am- 
plitude of the static paraniagnetic susceptibility of EuS, 
Co = (gd/Kp(0))2 = 0.19, Efgd = 0.24(2) was ob- 
tained. 

We also note that is smaller by a factor of four than 
Kg following from the longitudinal susceptibility and the 
measured value of Cq. We relate this difference to the 
dipolar interaction, which invalidates the proportionality 
^ E!~2{q), as can be inferred from the results of 
Refs.pUigf 



C. Linewidths 

The linewidths evaluated here are deflned by the 
Lorentzian shape which we assumed when analyzing the 
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spectra by Eq.(||L Note that already previous studies 
on powdered EuSEl favored this shape over the Gaussian 
and the truncatjed Lorentzian forms. Moreover, also the 
MMC approachll3 determined the damping of the mag- 
netization modes investigated here by assuming an ex- 
ponential relaxation at long times, which corresponds to 
the Lorentzian shape at not too high frequencies. 

The results for the widths of the spin wave peaks and of 
the central peak of the longitudinal fluctuations are dis- 
played by Fig.^ Within the experimental errors, there 
is no temperature variation down to the lowest temper- 
atures in region DC (see shaded area in Fig.^). Note 
that the errors for F^ are larger than those for Tgw be- 
cause of the smaller spectral weight of the longitudinal 
fluctuations, as it is seen in Fig.^. As the most striking 
result we infer from the presentations against q^{< q'^) 
in Fig.|^, that (i) the relaxation rates of both modes obey 
the simple relations 

r,„(g,T) = L,,„(^)2+r,^(0,T), (8a) 

<ld 

T,{q,T) = L,(^)2-)-r,(o,r), (8b) 

being indicated as solid lines in Fig.|^, and (ii) that the 
resulting kinetic coefficients Lsw = 35(2) ^eV and = 
40(3) fjieV agree within their uncertainties. Moreover, 
the damping of both modes remains finite with small val- 
ues of this background damping, Tsw{o^T) = 0.6(4) ^eV 
and Tz{q = 0) = 2.8(4) ^eV , which will also be discussed 
in the next section. 
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V. DISCUSSION 

We start from a rather general aspect of the critical 
phenomena, i.e. the scaling hypothesis extended to dy- 
namical quantities, like the relaxation rate of the order 
parameterB. Then the dipolar interaction in Heisenberg 
ferromagnets can fully be taken into accouat by intro- 
ducing q/qd as a second scaling variablellil'tj in the ho- 
mogeneous scaling function for the linewidths of the 
transverse modes investigated here, see Fig.|]b: 



Tl{q,T)^q^-g^ 



(9) 



q_ q_ 

Let us first approach the problem from the •param- 
agnetic side to which we designate the index \x — 
There according to both experimentEl and MMC theorytj 
the crossover from exchange dominated dynamics, being 
characterized by the exponent z = (Z) + 2)/2 = 5/2, to 
the dipolar critical dynamics with z — 1 occurs deep in 
the DC region, i.e. for g, k(T) ^ g^- This is the reason, 
why dipjaj.a|t,effects were not realized in the early investi- 
gationaiHa'cIa Once the dipolar dynamics, i.e. z = 2, 
has taken over very close to Tp, where Kp(T) ^ qd, 
the paramagnetic scalin g fu nction assumes the constant 
value (see. Eq.(7) of Ref|l|): 

5p(cx),0) = Ldg^2^ (10) 



FIG. 7: q-dependence of the Lorentzian linewidths of the a) 
spin waves and b) longitudinal fluctuations at selected tem- 
peratures measured in the dipolar critical regime, see Fig.|l^. 
For comparison are indicated data taken on powdered EuS at 
T = Tc. 



For EuS, the (dipolar) kinetic coefficient Ld has been 
determined from the analysis of the relaxation rates p£ 
the transverse-fluctuations in the limit q = abovecEl 
and below TjB: 



rj(g 



0,T) = 



(11) 



The value hLd = 38(4) fieV was iopjxL to agree 
fairly well with the MMC estimatetj'oEJ HLd = 

h.lgL^iB^kBTc|^lQ qPl^T^'' - 30 y^eV. 

Below Tc, as a matter of fact, we recover here the 
same g^-dependence of both relaxation rates F^ and T sw , 
Eq.(H), and in particular, within the experimental errors, 
their coefficients agree with Ld-, Lsw ~ Lz ^ Ld- These 
are the central results of our work. In terms of the dy- 
namical scaling hypothesis they imply that by passing 
the Curie temperature from above, the dipolar dynamic 
universality class for the transverse (critical) fluctuations 
characterized by z = 2 is not changed. This basic 
feature is nicely confirmed by (rather old) Fj(q) data 
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measured slightly above Tc of powdered EuS, which we 
have added to Fig.^a. Within the given error bars, they 
agree with respect to magnitude and g^-variation with 
the present data taken below Tc- We also observe 
a slight systematic enhancement of r^{q,T < Tc) over 
rj(g, Tc), which seems to be associated with the (small) 
background Tz^tg to be discussed below. 

To date, the fundamental phenomenon of dynamical 
scaling in Heisenberg ferromagnets has been established 
for the exchange ('true')critical regime (EC in Figjl^), 
i.e. for q > k{T) > qd, where F^ becomes independent 
of temperature, l|e.p|J|II] ^ (q/Qci)^^^ is maintained on 
both sides of TcaufiQiifl. Our experiment provides ev- 
idence for the 'true' (temperature independent) dipolar 
critical behavior extending to even for smaller q-values, 
q < Kz{T), see shaded region in Figjl^. A tempera- 
ture variation might set in when the dipolar hydrody- 
namic (DH) is reached. This extended dynamic critical 
behavior with z = 2 below Tc is also very much differ- 
ent from the situation above Tp-. where z = 2 is attained 
for extremely small q ^ qd QcJ, and only a small part 
of the static DC region is occupied by dynamic dipolar 
criticality. Both features are consistent with the general 
fact that the dynamic critical behavior reflects rmre de- 
tails of the system, in particular conservation lawsO, than 
the static properties, like the susceptibility tensor x(q), 
which define the critical regimes displayed in Fig.|l|a. 

Having established the dynamic class, z = 2, in a large 
part of DC below Tc, we are now able to discuss the 
further consequences of the q- variation of the damping, 
Eq. (H) . The kinetic coefficients and Lsw turned out to 
be identical for both the spin wave and the longitudinal 
modes. By looking at Fig.|lb, this result emerges from the 
continuity of the critical behavior of the transverse modes 
(5S^i(q) and (5S^2(q). Note that above Tc, the designa- 
tions 'transverse' and 'longitudinal' define orientations of 
(5S with respect to q being indicated by a superscript in 
(5S". Below Tc, due to the symmetry breaking through 
MsBz, this 'conventional' definition of the mode polar- 
ization (5S(q) is possible only for certain q-directions, 
like the one chosen in the present work. The 'trans- 
verse' modes drive the ferromagnetic transition on the 
paramagnetic side, where the q- vector clearly determines 
the symmetry. By passing Tc from above, Fig.nb sug- 
gests the continuous transformations, SS^^ 
(5S-^^ SS^ . Note that this is only true for our exper- 
imental configuration, where the order parameter Mg is 
oriented perpendicular to q. The fact that the kinetic co- 
efficients below Tc, Lsw and Lz are equal and, nw^over, 



agree with Lp = 38(2)/xeT^ of 5S (q) above Tc 



(see 



also Fig.0, can be immediately be related to the modes 
with transverse polarization 58 _L q, which arc the only 
ones to become critical in real Heisenberg ferromagnets, 
i.e. with dipolar interaction. On the paramagnetic side, 
these critical modes display relaxational dynamics, but 
below Tc their dynamical shape depends on the direc- 
tion of their propagation vector with respect to the order 
parameter M^. If these transverse critical modes propa- 



gate, for example, along Mg tii^ are predicted to exhibit 
spin wave dynamics in DCli3'E3 with pure pGoldstone- 
like susceptibilities X^^iSi II ^s) = ild/q)'^ EJ. The re- 
maining longitudinal mode, JS^ || q, should be strongly 
damped and suppressed in intensity. These modes can 
be studied in a configuration, where the polarizing field 
B is oriented parallel to the scattering plane. 

The dipolar symmetry with respect to q is also re- 
flected by the fact that in DC the direction of q is paral- 
lel to-the largest eigenvector of the susceptibility tensor 
x(q)ll2l, V3, while the second largest, V2, is parallel to 
Ms. This symmetry of x(q) changes if the DC-regime 
is left. Then the directions of these two eigenvectors are 
just interchanged for our configuration q _L M^, while 
for a general orientation between q and a gradual 
rotation, V3 — > V2 and V2 — > — V3, takes place. Hence, 
outside of DC, the 'leading' static symmetry of the fer- 
romagnet is defined by the order parameter Ms and, 
moreover, the dynamical exponent attains its isotropic 
value, z ~ 5/2. This implies that there both spin wave 
modes, SS'^^ _L q and SS^^ \\ q, are critical. Re- 
cently, their frequencies have-ibeen investigated in some 
detail by polarized neutronsEll. For the damping similar 
data are stilLlacking. Early work employing unpolar- 
ized neutronsci away from Tc of EuO provided relaxation 

rates, which where consistent with Tswiq,T)s^ ^^'^q^- 
This agrees with the scaling hypothesis, Eq.(g), provided 
gsw — [q/ 'iz{T))'^^ '^Osw {0,00). Unlike our observation in 
the DC region, in the exchange critical regime the dom- 
inance of the thermally excited spin waves gives rise to 
the strong increase of their linewidths with temperature. 
The leading q^-dependence has been predicted by Vaks 
et al.Lj and results from spin wave - spin wave scattering. 

Another interesting point is the fact, that both 
linewidths do not change with temperature down to the 
lowest values studied here, T/Tc — 0.92. This is some- 
what surprising wi| th | r egard to the existing results of the 
MMC-calculationsll§E3, which using some interpolation 
have also been indicated in Fig.^. We notice that they 
predict a slight temperature variation, i.e. a narrowing 
with decreasing T which is not observed. Also the q^- 
variations of our linewidths are not reproduced and the 
absolute MMC values for Fs^, and F^ are smaller and 
larger, respectively. This seems to indicate that the as- 
sumptions of the MMC approach, like the Lorentzian ap- 
proximation for the modes with large q and some cutoff of 
the dynamics, may not be valid. By a more phenomeno- 
logical point of view, we rather conjecture that the non- 
critical longitudinal fiuctuations, 6Sz \\ q, play a much 
more important role than believed to date. Their damp- 
ing is expected to be F^(q) = Ld/Xz — ^d, i-e. temper- 
ature independent and rather large and may provide an 
efficient relaxation channel for the critical modes. The 
importance of F^ has been realized recently in the relax- 
atioe. rate of the homogeneous magnetization F^ below 
Tc^. 

A final comment on the 'background' damping of the 
longitudinal fiuctuations, Tz{q = 0,T < Tc) = 2.8^ieV 
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may be appropriate. We conjecture here that - simi- 
lar as above Ea - damping effects by the longitudi- 
nal polarization of the SS^ mode are picked up within 
the resolution of our experiment. Due to the dipo- 
lar demagnetization, tliese modes are uncritical ( dy- 
namic exponent z=0 EJ). This conjecture is based on 
the fact that our background value is rather close to 
r.,bg = Lbg/xHl = 0,T < Tc) = 1.8{2)iieV, where 
we have inserted (i) the background kinetic coefficient 
Lhg = 1.8{2)fieV detecmined from the critical behavior 
of T^{q = 0, T > Tc jB, and (h) the susceptibility of the 
longitudinal spin wave modes, Xswi^ 0,T < Tc) — 1, 
which at T = Tq transform into the longitudinal param- 
agnetic modes x^il ^0). It may be interesting to note 
that this background seems to be absent or at least sig- 
nificantly be reduced in the paramagnetic relaxation rate 
rj(q, T) of Ref.|, depicted in Fig.0^ This indicates that 
the appearance of the order parameter increases r^(0, T) 
while its effect on the background rate of the spin wave 
linewidth, r^„(0,T) = 0.5(4) ^eV turns out to be small. 

VI. SUMMARY AND CONCLUSIONS 

We have conducted a study of the dynamics of spin 
fluctuations in the ferromagnetic state of EuS close to 
Tc, where dipolar effects are expected to play a domi- 
nant role. For intensity reasons we concentrated on the 
spin wave and longitudinal fluctuations, (5Sj^ and (5S^, 
both transverse with respect to the momentum transfer 
q, which was chosen perpendicular to M. This configura- 
tion allowed to extract from the spin wave frequencies the 
characteristic dipolar wavenumbcr qd and from the static 
susceptibility, Xz{^.^T) the correlation length — 
of the longitudinal fluctuations. The latter exhibits devi- 
ations from the Ornstein-Zernicke law at g <C Kz, which 
appear to be related to the pr|edicted thermal renormal- 
ization of Xz by the spin waveal3o. This feature, though 
subject of rather intense research in the past years (see 
e.g. References [35|j3^ , ^ , has not been identifled before. 
The reason for the observation of the mass renormaliza- 
tion of Xz (q, T) on EuS can be attributed to the rather 
large critical amplitude Kz{0) = 0.9lA~^, which allows to 
explore the regime q/Kz ^ 1- This is in contrast to the 



situation in the itinerant ferromagnet NiE3, where Kz(0) 
is almost one order of magnitude smaller. 

As the central result of our study, we regard the 
linewidths of both modes. They display the same q- 
variation — L^^q/ qdf' , which we explained by us- 
ing the dynamical scaling hypothesis and existing data 
for F^(q,T > Tq) of EuS. Apart from a small flnite 
background for Tz{q — 0), which was already observed 
above Tc in a previous work, the absolute values of both 
linewidths prove to be identical. Moreover, the relevant 
kinetic coefficient Ld agrees with the value obtained from 
the linewidths of the transverse fluctuations (5S J (q) mea- 
sured in the dipolar critical regime slightly above Tc- 
This quantitative feature suggests the transformations 
(5Spi SS'-^^, and (58^^ ^ (5S^ when passing Tc from the 
para- to the ferromagnetic side without any change of the 
dynamics and thus obeying dynamical scaling - as conjec- 
tured in Ref.|2^. As a surprising feature we note that the 
'true', i.e. temperature independent dipolar critical be- 
havior extends so far into the static DC region explored 
here, see Figj^a. This contrasts to the extremely narrow 
regime just above Tc, in which the dipolar anisotropy 
changes the dynamic universality class of the transverse 
modes from z = 5/2 to z = 2. The observed q- and tem- 
perature variations of our linewidths are notifuUy consis- 
tent with predictions by the MMC-theoryll3, which is a 
bit unexpected regardipgr-the success of this approach on 
the paramagnetic side.tilll3 

Finally we should note that for intensity reasons we 
were not yet able to measure the crossover of FJ to the 
hydrodynamic {q < Kz/9) and to the exchange critical 
(g > qd) regime, fo r whi i ch , sam e more theoretical work 
has been published. tilc3E3'HE3 An even greater exper- 
imental challenge is th|e-|low intensity of the longitudi- 
nal polarizations i5S^^E3 and 5S^. Their non-critical 
dynamics may be responsible for the absence of any 
temperature effects on the both linewidths investigated 
here. Other unsettled problems are the mechanisms re- 
sponsible for the small and non-critical background ef- 
fects occuring in the static susceptibilities, spin wave fre- 
quencies and relaxation rates. The values quoted here 
for Kg, Ksn],^z{0,T) and Fj^(0,r) may help to answer 
the question whether there exists a common origin, like 
anisotropy or flnite internal magnetic fleld. 
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